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The Closed Cluster method (CC method) is applied to find solutions for various calculation 
problems of the energy band structure of graphene. The essence of the CC method consists in the 
addition of closing bonds between edge atoms to the usual cluster method in order to eliminate the 
’’dangling” bonds on the edges of the cluster. We study the cases of an ’’infinite” layer of graphene 
as well as nanoribbons, nanotubes and bilayer graphene. Results for these cases are in agreement 
to that what was obtained by means of other methods (tight binding approximation and others). 
By means of the CC method we also study the problem of point defects in graphene and obtain 
the distortion of the energy spectrum. The energy spectrum of the layer Ci- X Si x (0 < x < 1) is 
found as well as the dependence of the energy gap on the concentration of silicon. We show that 
the energy band structure of Ci- X Si* looks like a tunnel transition. Wave functions of graphene in 
the symmetry points of Brillouin zone are also obtained. 
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I. INTRODUCTION 

Studies of various properties of graphene and it‘s appli¬ 
cations have attracted much attention in the last years, 
as it is well-known [1]. In this paper, it is proposed for 
an approach based on the use of closed clusters (CC) 
to calculate the energy band structure of graphene. We 
have developed this approach earlier in application to the 
one-dimensional and three-dimensional crystals with di¬ 
amond structures [2]. The essence of the CC approach 
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is to bring together all bonds of atoms which are located 
on the edge of a cluster in order to eliminate the “dan¬ 
gling” bonds. This approach is found to be especially 
useful for calculation of the energy spectrum of crystals 
with point defects, such as vacancies or impurity atoms. 
The simplest and most widespread approach to calculate 
such impurity states is known to be the effective mass 
method. This approximation works sufficiently well in 
cases of impurity levels being located along the borders 
of the energy zones. However, it is inapplicable to the 
description of so-called “deep levels”, which lie far from 
the zone borders [3]. One of the methods used for those 
levels is the cluster approach, in which a group of atoms 
- the cluster - is mentally picked from a crystal lattice. 

This cluster is considered as a separate “molecule” and 
for it‘s calculation the usual quantum chemistry methods 
are applied. The advantage of the clusters approach is 
it's applicability for modelling real situations of impurity 
atoms and - if necessary - taking into consideration a 
possible distortion of the crystal lattice. 

However, the usual cluster approach has a deficiency. 
Is the “infinity” crystal substituted by a group of atoms, 
a problem with atoms lying on the edge of the cluster 
arises. The presence of such atoms with torn bonds dis¬ 
torts the energy spectrum of the crystal. This distortion 
can be diminished by increasing the size of the cluster, 
but the approximation to the exact value is very slow. 

This deficiency can instead be removed by connecting 
the torn bonds with each other and hereby closing them. 
A similar procedure, the so-called “periodical boundary 
conditions”, is applied in studying the energy spectrum 
of the infinite ideal crystals. The special feature of our 
approach is the application the closing procedure of these 
bonds to small clusters to study crystals with distorted 
regularity properties (such as impurity, edges and other). 
In this paper, the CC approach is applied to graphene, 
a relatively new material with numerous of unique prop¬ 
erties [4]. In section II the fundamental idea of the CC 
approach is presented as well as the rules for the construc¬ 
tion of diagrams corresponding to the various clusters. 
In section III examples of building Hamiltonian matrices 
and solutions for the secular equations in case of period- 
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ical structures without edges and defects are given. The 
aim of this section is to test our approach by means of 
comparison of the results with those obtained by other 
methods. In section IV the CC approach is applied to 
bilayer graphene. In section V we present wave functions 
at symmetry points of the Brillouin zone. In section VI 
we study nanoribbons and nanotubes by means of closed 
clusters. Section VII contains the calculations of the 
graphene energy spectrum in the presence of impurity. 
Finally, in VIII we study the energy band structure of 
the hypothetic monoatomic layer Ci-^ Si x (0 < x < 1). 
Section IX concludes with a discussion of the obtained re¬ 
sults and furthermore describes a possibility for making 
the approach more precise as well as the application of 
the CC approach to other problems. 

II. THE ESSENCE OF THE CC APPROACH 

AND A RULE FOR CONSTRUCTING THE 
DIAGRAMS 

By the term “closed cluster” we mean a group of N 
atoms which reflects the structure of the crystal and fur¬ 
thermore the atom bonds which are present within the 
group - at the same time the torn bonds are absent. A 
similar procedure of closing the bonds known as “periodi¬ 
cal boundary conditions’ is often used in one-dimensional 
atom chains when atom N+l is identical to the first one. 
Such a procedure is also possible for 3D-systems, but is 
rarely applied here, since closing a 3D cluster is much 
more difficult. In 2D cases however, closing can be re¬ 
alized very easily - therefore, application of the CC ap¬ 
proach to graphene and other 2D structures seems to be 
very effective. 

While constructing clusters corresponding to graphene 
one must start out from it‘s crystal structure. The latter 
is well-known, a hexagonal layer which is formed by two 
sublattices, here A and B [5]. The simplest unclosed 
cluster corresponding to such a structure is shown in 
Fig. 1(a). 

The main idea of the closed cluster approach is the 
necessity of closing the torn bonds shown in Fig. 1(a) in 
the way shown in Fig. 1(b). In case of an endless graphene 
layer the structure of the lattice requires the following 
rules regarding the construction of a closed cluster with 
approximation of the nearest neighbors. 

• Rule 1. Each atom of sublattice A is bonded with 
three atoms of sublattice B and vice versa. 

• Rule 2. All real and closing bonds have the same 
energy levels. 

• Rule 3. The number of atoms in a cluster N must 
be divisible by six. 

The last rule is explained in detail in section V. 

All three rules in fact are satisfied within the construc¬ 
tion in Fig. 1(b), although the second rule seems to be 
broken at first sight: for example, the closing bond be¬ 
tween atoms 1 and 4 looks different from that between 1-2 



FIG. 1: The most simple clusters of the graphene lattice 
(a) unclosed cluster N = 6 
(b) closed cluster N = 6 
(c) cluster surrounded by neighbor atoms 
dashed line: closing bonds 
continuous line: real bonds in graphene lattice 


or 1-6, since atom 4 is further away from 1. However, it 
must be kept in mind that the bond 1-4 is not a real bond. 
In distinction from Fig. 1(c), the closed cluster 1(b) must 
be understood as a diagram or graph, which makes it 
easier to obtain the Hamiltonian matrix. Furthermore, if 
these Hamiltonian matrix elements corresponding to the 
closing bonds are chosen the same as for existing bonds, 
then the “interaction” between atoms like 1 and 4 in 1(b) 
in fact describes the interaction between 1 and 1’ in 1(c), 
latter of which is absent in cluster N = 6. 

It must be noted, that for the numbers of atoms in 
sublattices A and B it is necessary to be equal, as closing 
bonds is only possible between atoms of different sublat¬ 
tices. For example, in cluster 1(b) it is only possible to 
close the bonds of atom 1 with atom 4, but not with 3 
or 5. This restriction is necessary to satisfy Rule 1. 

We also note that energy values obtained from the 
solution of secular equations are independent from the 
choice of the cluster in case of an “infinite” ideal lattice, 
but are only dependent from the number N. They are 
as well independent from the numeration of the atoms, 
since changes in numeration only cause determinant per¬ 
mutations. 

Fig. 2 shows some of the possible clusters with N = 12 
with dissimilar ways of closing. 

To conclude this section, we note, that the name 
’’closed cluster” which we use is different from the term 
’’closed walk” used in graph theory [6]. For example, the 
cluster shown in Fig. 1(a) is a ’’closed walk”, but not a 
closed cluster - the latter is shown in Fig. 1(b). 

One definition of a closed cluster can be given as a 
cluster, where each atom is linked to the same number of 











1, if atom m is bound with atom n 
0, if atom m is not bound with atom n, 
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FIG. 2: Examples of closed clusters for N = 12. 


neighbor atoms as in the corresponding crystal. 


V mn — 


E-E 0 

e =- 

7o 


with 

E - energy of an electron 

Eq - energy of an electron in the |2 p z ) state in an iso¬ 
lated carbon atom. We use Eq = 0 

and 

—70 - the transfer integral between neighbor atoms 
with wave functions with the same sign, 70 > 0. 

The problem of calculating the corresponding values 
of energy 6 (in 70 units) is then reduced to solving the 
secular equation 


det M mn = 0. (4) 

Now we consider some examples of applications of the 
approach. We begin with the simplest closed cluster 
N = 6, which is shown in Fig. 1(b). The secular equation 
(4) has the form (empty matrix cells standing here and 
further for zeroes): 


III. BUILDING OF A HAMILTONIAN MATRIX 
AND SOLUTION OF SECULAR EQUATIONS 

The basic idea underlying the CC approach is the same 
as in the usual molecular orbital approach (MO) [7]. In 
particular, this is the representation of the wave function 
of the cluster T(r) as a linear combination of the wave 
functions of the atoms 


N 

y(r) = J2c n <p(r-a n ). ( 1 ) 

n= 1 


Functions p in graphene are | p z ) - orbitals of carbon 
atoms, with axis z being perpendicular to the layer. We 
designate wave functions along and opposed this axis as 
”+” and 

The standard procedure of obtaining secular equations 
leads to a system of linear equations which help to find 
the coefficients c n in 


N 

^ ^ MmnCn = 0? 771 = 1, 2, . . . , A, 

n= 1 

where 


Mmn — mn T (1 


^mn^V 


mm 


( 2 ) 

( 3 ) 


(5 


mn 


1, 771 = 71 
0, 771 7^ 71, 


e 1 1 

1 5 1 1 


1 1 5 

1 1 1 

Calculating the determinant (5) leads to the following 
equation for finding 5: 

£ 4 (£ 2 -9)=0. (6) 

The solutions of (6) ±3,0,0,0,0 are exactly those en¬ 
ergy values at bottom, top and Dirac points of two en¬ 
ergy bands, which are obtained from a tight binding and 
nearest neighbor approximations [5]. 

The fact, that the CC approach provides exact values 
of energy band boundaries already at a minimal cluster 
size is very important in the calculation of the energy of 
impurity states, as the energy of such states are counted 
from the bands’ boundaries. 

Let us consider the next example with N = 12. If we 
proceed from the clusters from Fig. 2, then the secular 
equation has the form 


= 0. 


( 5 ) 

























4 


D 12 = 


£ 1 
1 £ 1 
1 £ 1 


1 £ 1 
1 £ 1 


1 


1 


1 £ 1 
1 £ 1 


1 £ 1 
1 £ 1 
1 £ 1 
1 £ 1 

1 1 £ 


= 0. 


( 7 ) 


The solution of equation (7) is given by following val¬ 
ues of energy e% ±3, ±2, ±2, ±1, 0, 0, 0, 0. 

The determinant of cluster N=24 has a form which is anal¬ 
ogous to (7) £>i 2 . The obtained energy values are shown in 
Table I. 

The main conclusion which may be drawn from the com¬ 
parison of cases N = 6 , 12 and 24 is the fact, that the energy 
values of the bottoms and tops of the lower (—3, 0) and upper 
band (0, 3) are equal in each case, therefore they are indepen¬ 
dent of the cluster size. With the growth of the cluster size 
new energy levels arise, however, previous levels remain the 
same. 

Each energy value in Tab. I is corresponding with some 
point of the Brillouin zone. More in-detail discussion on that 
issue is presented in sec. VAt this point we just briefly discuss 
the fourfold degeneration of level £ = 0 at all N. 

At first glance it seems not to be in accordance with £ being 
0 at six corners of the Brillouin zone (points Ah, K 2 , •.., Kq 
in Fig. 3). 



FIG. 3: Brillouin zone of graphene [5]. 


section we use a three-dimensional diagram to study bilayer 
graphene. 


IV. BILAYER GRAPHENE 


For N = 12, the 3D closed cluster corresponding to bilayer 
graphene is shown in Fig. 4. In addition to the closing bonds 
within the layers (energy 70 ), also the bonds describing the 
interactions between those layers (energy 71 in 70 -units) are 
present in the cluster. 

Note, that Fig. 4 does not show the displacement of the 
upper layer towards the lower, since this Figure doesnt repro¬ 
duce a real atomic structure, but is solely a diagram which is 
used for building a Hamiltonian matrix. 



- real bonds within the layers y 0 , 

. closing bonds within the layers y 0 , 

... bonds between layers yj. 

FIG. 4: Cluster N = 12 for bilayer graphene. 

The secular equation corresponding to Fig. 4 has the form 


D 


( 2 ) _ 


1 


71 


1 


£ 1 
1 £ 1 
1 £ 1 


71 


71 


71 


71 


71 


71 


1 £ 1 
1 £ 1 


71 


71 


1 £ 71 

71 ^ 1 

1 £ 1 


1 £ 1 
1 £ 1 


71 


1 £ 1 
1 £ 


= 0. 


Solving equation ( 8 ) leads to following energy levels: 


(8) 


However, it must be kept in mind, that only two of these 
six points, K\ and K 2 for example, belong to the first Bril¬ 
louin zone, others belonging to the next zones. The twofold 
degeneration at points K\ and K 2 is what leads to fourfold 
degeneration of the level £ = 0 . 

In conclusion of this section it should be noticed, that the 
rules for diagram construction formulated in sec. Ill does not 
require these diagrams to be plane. Therefore, they also 
can be applied to nanotubes as well as spherical surfaces 
(fullerene). Stratified cluster can also be used. For example, 
the diagram in Fig. 2(a) can be rearranged to a hexahedron 
without breaking the bonds. The application for the space di¬ 
agram is useful for studying many-layer graphene. In the next 


£ = ±(3 + 71 ), ±(3 - 71 ), ± 71 , ± 71 , ± 71 , ± 71 . (9) 

The main result that ensues from (9) is that there is an 
energy gap e g — 27 . The presence of this gap in bilayer 
graphene is a fact discussed in many papers [ 8 ] [9]. In the 
model which we are using, the energy gap is caused by the 
interaction between the layers and disappears if we set 71 = 0 . 


V. WAVE FUNCTIONS 

Within our considered model the cluster wave functions are 
determined by the totality of coefficients c n , according to ( 1 ). 
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TABLE I: Electron energy values with varying numbers of atoms in the cluster 


N 

Energy 

6 

± 3 










0 

0 

0 

0 

12 

± 3 



± 2 

± 2 


± 1 




0 

0 

0 

0 

24 

± 3 

±(\/3+l) 

±(\/3+l) 

± 2 

± 2 

± 1 

± 1 

± 1 

±(\/3-l) 

±(\/3-l) 

0 

0 

0 

0 


To obtain the latters one must solve the system of equations 
(2) under defined energy values e. For £ = 3 and e = — 3 the 
result is obvious and shown in Tab. II. Let us now consider 
the cases e — ±1 and e — 0 in more detail. Latter is most 
interesting, since at this point the valence band meets the 
conduction band. 

As an example, let us consider the cluster N — 12 , which 
is shown in Fig. 2(b), but with a different atoms numbering 
(Fig. 5). 


• * 


Cl - C 7 

C2 = Cq = C 8 = Ci 2 
C3 — C5 — Cg — Cn 
C4 — Cio- 


( 11 ) 


Then the task is to solve just four equations instead of 
twelve. 

If we take (11), the equations have the form 



FIG. 5: Cluster N = 12 for calculation of graphene 
wave functions. 


To determine the coefficients ci, C 2 ,..., C 12 in (2) we must 
solve a system of twelve linear equations. For cluster in Fig.5 
it has the form 


(1) 

6 Ci + C 2 + Cq 

+ 

C 7 

0 

(2) 

Cl + £C 2 + C 3 

+ 

eg — 

0 

( 3 ) 

C2 + SC 3 + C4 

+ 

eg = 

0 

( 4 ) 

C 3 + £C 4 + C 5 

+ 

cio = 

= 0 

( 5 ) 

C4 + £C5 + Cq 

+ 

C11 = 

= 0 

( 6 ) 

CI + C 5 + £Cq 

+ 

C12 = 

= 0 

( 7 ) 

Ci + £C 7 + Cg 

+ 

C12 = 

= 0 

( 8 ) 

C 2 + C 7 + £C 8 

+ 

Cg — 

0 

( 9 ) 

C3 + Cg + £Cg 

+ 

Cio = 

= 0 

(10) 

C4 + Cg + £Cio 

+ cn 

= 0 

(11) 

C 5 + Cio + SC\\ 

. + C12 = 0 

(12) 

Cq T C7 T Ci 

1 + SC12 

= 0 


As we can see from Fig. 5 our cluster is reflection-symmetric 
in the axes y and x, due to the closing bonds 1-6 and 7- 
12 . Hence, wave functions must be symmetric (S x ,S y ) or 
antisymmetric (A x A y ) to reflection in the x and y axis. All 
together there are four possible symmetries of wave functions: 
S x S y , S x A y , A x S y and A x A y . 

First we consider the S x S y case. Following relations among 
the coefficients must be satisfied 


l.(e + l)ci + 2 c 2 

= 0 

2 . 

Cl + (5 + 1 )C 2 + C3 

= 0 

3. 

C 2 + (s + 1)C3 + C4 

ss= 0 

4. 

2 c 3 + (c + l)c4 

= 0. 


( 12 ) 


The determinant of (12) has to be zero and therefore we 
obtain 


S x S y :e = -3,-2,0,l. (13) 

In other symmetry cases we have 


S X A 


y 


A x Sy 


A x A y 


ci — — C 7 , 

C2 — Cq — —C 8 = —C 12 , 
C 3 = C 5 = — Cg — — Cn, 
C 4 = —C 10 . 

Cl = C 7 = 0, 

C 2 = —Cq m Cg = —C 12 , 
C 3 = —C 5 — Cq — — Cn, 
C 4 = Cio = 0 . 

Ci — c 7 — 0 , 

C2 = —Cq = — Cg — C12, 
C3 = —C5 — —Cg — Cn, 
C4 = Cio = 0 . 


(14) 


(15) 


(16) 


In the latter two cases there is ci = C 4 = C 7 = cio — 0 
due to the antisymmetry to axis x (A x ). Therefore, only two 
equations remain in each case, namely those of the coefficients 
C 2 and C 3 , and consequently only two values of energy for each 
A x S y and A x A y . 

The energy levels for various symmetries are shown 
in Tab. II, Tab. Ill shows the values of the coefficients 
ci, C 2 ,..., C 12 . These coefficient values are obtained by solv¬ 
ing the system of equations like ( 12 ) and analogous systems 
for other symmetries for the energies e = 0. ± 1, ±3. The set 
of these coefficients determines the wave functions for various 
symmetries within our considered model. Tab. II and III also 
show the Brillouin zone symmetry points, which can also be 
seen in Fig. 3. 

For the purpose of illustration one can represent these func¬ 
tions graphically by extending the cluster over the whole 
graphene layer. In the case e — ±1 the functions coincide 
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TABLE II: Energy levels of cluster N = 12 corresponding to various symmetries of wave functions 



TABLE III: Wave functions in the symmetry points of the Brillouin zone 


^Energy 

Cn 

-3 

-1 

0 

0 

0 

0 

1 

3 

Cl 

1 

1 

1 

1 

0 

0 

1 

1 

C 2 

1 

1 

- V 2 

- 1 / 2 

1 

1 

-1 

-1 

c 3 

1 

1 

- V 2 

- 1 / 2 

-1 

1 

1 

1 

c 4 

1 

1 

1 

-1 

0 

0 

-1 

-1 

C5 

1 

1 

- 1 / 2 

- 1 / 2 

1 

-1 

1 

1 

c 6 

1 

1 

- 1 / 2 

V 2 

-1 

-1 

-1 

-1 

C7 

1 

-1 

1 

-1 

0 

0 

1 

-1 

C8 

1 

-1 

- 1 / 2 

V 2 

1 

-1 

-1 

1 

c 9 

1 

-1 

- 1 / 2 

V 2 

-1 

-1 

1 

-1 

ClO 

1 

-1 

1 

1 

0 

0 

-1 

1 

cn 

1 

-1 

- 1 / 2 

V 2 

1 

1 

1 

-1 

C 12 

1 

-1 

- 1 / 2 

V 2 

-1 

1 

-1 

1 

Symmetry 

Point of BZ 

S X Sy 

r 

SxAy 

M - 0 

S X Sy 

K\ 

S X Ay 

k 2 

A x Sy 

K] 

Ax Ay 

k 2 

SxSy 
M + 0 

S X Ay 

D 




y y 




FIG. 6: Wave functions of graphene various symmetries 
at the points K. The data is taken from Table III. 


with those in [5]. For the case e = 0, wave functions for 
various symmetries are shown in Fig. 6. 

From Fig. 6 it is easy to see, that in cases 6(a) and 6(b) 
the wave functions are plane waves with A = | a, where a is 
the length of the vector connecting the nearest neighbors of 
atoms A and B. The fronts of these waves are parallel to axis 
x , which means that the wave vector k is directed along the 
axis y. 

The length of the wave vector of these waves k = ~ 
coincides with the length of k in Brillouin zone point K\ 
(Fig. 3) - thus, Figures 6(a) and 6(b) show the wave functions 
at point K\. In cases 6(c) and 6(d) one can easily see, that 
the fronts of the wave of each sublattice are at an 120° angle 
to the x axis. Therefore, the waves vectors are directed at a 
30° angle to the x axis, that is the point K 2 in Fig. 3. The 
length of these waves is also § a. As a consequence, Fig 6(c) 
and Fig. 6(d) show the wave functions at point K 2 . 

The fact that all four wave functions shown in Fig. 6 cor¬ 
respond to the energy e — 0 can be directly deduced from 
the form of these functions. Namely, the total energy of the 
interactions of the atoms with their respective nearest neigh¬ 
bors is equal to zero, taking into consideration the values and 
signs of c n from Fig. 6. 

In Fig. 6 it is also apparent, that those structure elements, 
which are periodically repeated within the wave functions at 
e — 0 are triplets of atoms: 1,—1,0 (Fig. 6(b)), 1, — T ? — T 
(Fig. 6(a)) etc. Consequently, the energy value e — 0 is only 
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(a) 



(b) Related 
closed cluster 
N=12 


FIG. 7: Armchair ribbon (a) and the corresponding 
closed cluster TV = 12 (b). 


possible in clusters, in which the number of atoms N is divis¬ 
ible by three. Since N also has to be even, due to the equal 
number of atoms in sublattice A and B, the number N finally 
has to be divisible by six. This is the base for rule 3 of cluster 
building from section II. 

It is easy to see that in clusters N = 8,10 and other N 
not divisible to six, energy values £ = 0 are absent. There¬ 
fore, only clusters with by six divisible N must be taken into 
consideration in the case of endless layers. 


VI. NANORIBBONS AND NANOTUBES 


The CC approach is especially applicable to impurities, as 
well as nanotubes and nanoribbons. In the latter two cases a 
group of atoms lying across the tubes or ribbons is a natural 
cluster, in which the number of atoms depends on the diam¬ 
eter of the tube respectively the width of the ribbon (Fig. 7). 
Let us now consider the armchair ribbon with N = 12. As 
zigzag nanoribbons are not analogous to carbon nanotubes 
[12] the former will not be discussed in this work. The above 
mentioned cluster is shown in (Fig. 7(b)). This cluster differs 
from the cluster in Fig. 2(b) in a single, yet very significant 
aspect: the cluster in Fig. 7(b) has no closing bonds between 
atoms 1-6 and 7-12. This reflects the fact, that the latter 
atoms are situated on the edge and therefore are connected 
to the ribbon by means of two, not three bonds. 

As one can see, the cluster in Fig. 7(b) belongs to a well 
known ladder type. The determinant Djv of this cluster may 
be written in block form for any even N 


Dn 


Ln 1 (e) Tvi 

//Vi Ln x (s) 


(17) 


where 



£ 1 
1 £ 1 
1 £ 1 


TjVi (£) = 


>N i (18) 


1 £ 1 
1 £ 


and T/Vi being the unity matrix of order Ah. 

The determinant (17) can be brought into a quasi- 
triangular form. For this, we need to add row number 
(Ah + 1) to the first row, row (Ah + 2) to the second and so 
on until row Ah is added to row 2TVi. Then, the first column 
is subtracted from column (Ah + 1), second from (Ah + 2) 
and analogously column Ah from 2Ah. 

As result we obtain 


Dn — Ln x (£ + l)TiVi ~ 1)? (19) 

where L n (x) are well-known polynomials, expressions of 
which are given in [10] (for g n (x) there) for n < 20. It is also 
known that representation of L n (x) via the trigonometrical 
functions is [11] 


where x — 2cos0. 

Then, for obtaining the energy values Dn(s) = 0 is reduced 
to conditions, which need to satisfy at least one of two of the 
following equations: 


Ln 1 (s + 1) — 0 

(21) 

or 


Ln 1 (£ — 1) = 0. 

In view of (20), these equations lead to 

(22) 

sin[(Ni + 1)0+] = 0 

and 

(23) 

sin[( Ah + 1)0-] = 0, 

where the 0+ are determined from 

(24) 

2 cos6± = £ d= 1. 

The solutions of (23) and (24) are 

(25) 

_ n±7r 

9± ~ N 1 + 1’ 

(26) 


where n± = 1,2 ,... ,Ah. 

Then, from (25) we get the formula for the energy values of 
a closed cluster with 2Ah atoms, which describes the nanorib¬ 
bon 

£n ± = 2 ^1 ± L ( 27 ) 

Taking into account that n+ and n_ both may possess Ah 
values and due to the ±1 part we see from (27), that s n± 














possesses 2N± = N values, as this has to be the case for a 
cluster with N atoms. 

If the ribbon is transformed to an infinite layer (iVi —> oo), 
then we get the boundary energy spectrum values from (27), 
which are £max = 3 at n+ = 1 and Emin — —3 at n_ = N ±. 
These values also coincide with the result in section III. 

Using equation (27) one can easily find the values of Ah, 
under which the solutions e n ± — 0 exist, that is the ribbon is 
metallic. The condition for this is cos = \ or cos = 

— |, otherwise — f, n++i ~ T"’ n + ~ 2n. Denoting 

n- — M one can write the both latter equations as 

Ni = 3M — 1, (28) 

with M being an integer. 

The formula (28) is the same as has been obtained by means 
of the tight binding model in [13] and [14]. 

The wave functions of nanoribbons with metallic conduc¬ 
tivity can be found by using the symmetry like it is done for 
the infinite layer in section V. 

We consider as an example the case N± — 5. The closed 
clusters for the description of the ribbon are shown in figures 
8 (a) and 8(b). According to the graph theory terminology [6], 
these clusters are isomorphic to each other as well as to the 
ladder graph shown in Fig. 8(c). The energy spectra of these 
graphs are therefore the same. 



FIG. 8: Wave functions of clusters for description of rib¬ 
bon Ni = 5 (a,b), ladder cluster (c) and ribbon with 
symmetry A x S y at 5 = 0 (d). 


Likewise to the case of infinite layer one can also obtain 
the energy values of the cluster, which are corresponding to 
various symmetries. Results for N\ = 5 are shown in Tab. IV. 

If only the energy spectrum is required it can be directly 
obtained from equation (27). 

The cases for other N\ can be considered analogously. Fig¬ 
ure 9 shows the wave functions of the symmetry A x S y for the 
Ni = 8 ribbon at e = 0 . 

If we compare the ribbons shown in Figures 8 (d) and 9 we 
can see the meaning of the metallic conductivity condition 
(28): when increasing Ah, in order to keep the energies of 
each atom at zero, it is necessary to add two rows of atoms, 
with positive and negative \p z ). The latter must however be 
separated from the previous row by including a row of zeroes. 
Altogether, we need to add three rows, which is where the 
term 3 M in (28) derives from. From Figures 8 and 9 one can 
also see that the distance between the nodes ^ = |a and the 
wave front is parallel to axis x. This corresponds to the wave 
vector k at point K\ of the Brillouin zone. It is easy to see, 
that in case of the symmetry A x A y the energy value e — 0 is 
corresponding to the point Kz , similar to the infinite layer. 

Using Ah is not the only way to describe the criterion for 
metallic conductivity. It also can be written by the means of 



FIG. 9: Wave functions of symmetry A x S y for the 
ribbon Ah = 8 at 5 = 0. 


the width of the ribbon L. Figures 8 (d), 9 and equation (28) 
show that it is possible to obtain the connection between M, 
Ah and L as in Table V. 

The relation between L and M in Table V can be written 
as 


L met = 3M_ 2 , Mat 1,2,3... . (29) 

This formula differs from L = 3 M + 1 in [15], which is 
maybe due to the different definition of the ribbon width in 

[15]. 

Figure (10) shows the dependence between the energy gap 
£ g and the ribbon width. Here it is e = 2ei, where £1 is the 
nearest level to zero and is found via (27). 

L 



FIG. 10: Dependence between the energy gap e g and 
the ribbon width. 


The values of Ah and L, where e g — 0 are not marked, they 
can be obtained from equations (28) and (29). 

This described approach of the calculation of nanoribbons 
by the means of ladder clusters may also be applied to nan¬ 
otubes, with some modifications. The formal transition from 
the nanoribbon to the nanotube consists in closing a cluster 
of the type shown in Fig. 7(b), so that closing happens be¬ 
tween atoms 1-6 and 7-12 - this way we get the cluster from, 
Fig. 2(b). If we wrap it up into a cylinder form, we obtain a 
cluster with three hexagons. It can be easily verified that in 
clusters of this type the number of hexagons in the cylinder 
section is always an integer. 

However, a cluster of type as in Fig. 8 cannot be wrapped up 
in such a way, because the atoms of the upper and lower edge 
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TABLE IV: Energy spectrum of a closed cluster for description of a Ni = 5 ribbon 


Symmetry 

Energy 

S X Sy 

—v/3-1 


-1 




V3-1 




S X Ay 




— a/3 + 1 




1 


a/3 + 1 

A x Sy 


-2 



0 






A X Ay 






0 



2 



TABLE V: The connection between the ribbon width L 
(scaled by a), Ni and M in the case of metallic 
conductivity. 


M 

JVi 

jjmet 

1 

2 

V 2 

2 

5 

2 

3 

8 

7 2 

4 

11 

5 





are from the same sublattice, A or B - and nearest neighbors 
always have to be from different sublattices. 

Therefore, in order to wrap up a cluster into a cylinder one 
must add a row of atoms. These atoms need to belong to a 
different sublattice than the atoms at the edges of the initial 
cluster. In an Ni = 5 example this means that the number 
N c of atoms situated along the circle cylinder section must 
equal six. 

It may also be easily verified that for clusters of both types 
N c must be even (which is a consequence of the equal number 
of A and B atoms) and the number of hexagons Ns along 
the circle has to be an integer. Hereat, beginning with the 
smallest value N s = 2 it is N c = 2 N s . 

In case of metallic conductivity (e g = 0) as well as in the 
infinite layer (section V), it is necessary for N c to contain an 
integer number of atomic triplets +, —, 0, that is, N c needs to 
be not only even but also divisible by three. All in all N c has 
to be divisible by six, and therefore we can write the criterion 
of metallic conductivity of nanotubes as 

N™ et = 6M, M = 1,2,3.... (30) 

For all other values of N c there is e 7^ 0. For example, in 
case of the least values N s = 2, iV c = 4 one must consider the 
closed cluster N = 8. It is possible to calculate a cluster of 
type from Figure 2(b), but for N = 8. 

Solving equation (4) for this case gives us the following 
energy spectrum: 


VII. IMPURITY IN GRAPHENE 

The cluster approach is, as already mentioned, especially 
appropriate for studying the impurities in crystals. That is 
connected with the fact that the impurity influence on the 
crystal depends mainly on the interaction with atoms, which 
are located in the nearest environment of this impurity. Ex¬ 
actly those atoms together with the impurity are the ones 
which can be considered as a cluster with N atoms. The 
study of this cluster is the essence of the impurity problem in 
the present approach. 

In this paper we restrict the consideration to the simplest 
cluster N = 6. Let us assume, that in the first lattice point 
of this cluster atom C is substituted by another one. This 
atom only differs from C in terms of the \p z )-electron energy 
of an isolated atom, which is Eo + AEo instead of just Eq. It 
is also assumed, that this impurity atom is silicon, although 
is also might be germanium or any other atom with the same 
outer electron shell structure as C. Furthermore, we assume 
the transfer integral 70 equal between all neighbor atoms. In 
order to account for the difference in these parameters within 
the cluster approach, one simply has to replace the ones in the 
corresponding Hamiltonian matrix elements with parameters, 
which characterize the transfer integral between the impurity 
and the surrounding atoms. The issue of the influence of these 
parameters on the energy spectrum in 3D crystals has been 
investigated by us earlier by means of Green’s functions [16]. 
For the purpose of this paper, however, this problem is not 
of a fundamental meaning, in so far as the main parameter is 

AE 0 . 

An example of a closed cluster with one impurity atom Si 
is shown in Fig. 11. 



FIG. 11: Closed cluster N = 6 with one impurity atom 

Si. 


£ c = ±1,±1,±1,±3. (31) 

Hence, the value of the energy gap for N c = 4 is e c g — 2, 
which is the maximum value of £ g for nanotubes. With an 
increasing N c the value of e c g decreases. 

It should be noted, that e c in (31) and e g are scaled by 
7o : £ c = , where 7§ is the tube’s transfer integral (hopping 

energy). This latter value is different from the value of 70 for 
the infinite layer and is dependent from N c , which should be 
taken into the account when s c g in eV needs to be obtained. 
However, this dependence does not exert any influence on the 
criterion e g — 0. 


The secular equation for this cluster has the form 


E 6 = 


(e - A) 1 
1 £ 
1 

1 

1 

1 


1 

1 

1 

1 

£ 1 

1 £ 


= 0, (32) 
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where A — 

70 

This sixth degree equation is reduced to a cubic equation, 
which can be solved analytically. As result we get following 
solutions 


£1,2,3 — 0 . 


£ 4 , 5,6 - 3 


A + 2 \J A 2 + 27 cos ( ^ + n ~^~ 


(33) 


where 


A(A»-f) 

COS if = —--. 

(A 2 + 27) 2 

It is easy to see that under A = 0 there is 

£4 — 3, £5 = —3, £q — 0, 
that is we obtain the result (6). 


(34) 


The numerical calculations of the dependence of e from A, 
which may be obtained from (33) and (34) or, simpler, directly 
from the calculation of the determinant (32), are shown in 
Figure 12. 


VIII. THE ENERGY BANDS OF 
MONOATOMIC LAYER Ci_* SI* (0 < x < 1) 

By means of the CC approach one may as well calculate 
the energy spectrum of a hypothetical object, a monoatomic 
layer of the type Ci_* Si*. That such an object can be created 
results from the following reasoning: by obtaining by means 
of epitaxi, graphene can be formed on the surface of SiC. The 
monolayer of Si, which appear hereat as a buffer layer, has 
also a hexagonal structure [17]. That is why by means of 
epitaxi, perhaps not only graphene can be formed but also 
the monoatomic layer of the type Ci_* Si* (0 < x < 1), with 
a hexagonal structure similar to graphene. 

In the study of this problem we restrict ourselves to con¬ 
sider only the simplest structure N — 6. This is found to be 
sufficient for describing in outline the energy band structure 
of Ci—* Si*. 

We consider the cases separately when in a cluster of 6 
C-atoms one, two, three, four or five atoms are replaced by 
Si-atoms. This corresponds to the values x — re- 

spectivaly. The case of x — | has been already considered in 
the preceding section. It should be noted that under ” impu¬ 
rity” we mean the atoms which are present in a lesser number. 
That is at x — | the impurity atoms are those of Si and at 
x = | those of C. 

The clusters, corresponding to these x-values, can be of two 
types. Namely, by x = | two impurity Si-atoms can interact 
with each other (type II) or not interact (type I). The same 
can be obtained for C-atoms by x = |. First we study the 
clusters of type /, which are shown in Fig. 13. 



A 

FIG. 12: Dependence of energy levels of cluster N = 6 
on A. The level 5 = 0 is threefold degenerated. 



16 16 16 1 

(a) x = 1/3 (b) x = 1/2 (c) x=2/3 (d) x 


I IO. 13: Clusters of type / with V = 6 for Ci_ a 



SO 


The secular equations for clusters shown in Fig. 13 can eas¬ 
ily be estimated. For example, by x = | we obtain 


(£ - A) 1 
1 e 


1 

1 

(e - A) 1 
1 £ 
1 

1 


1 

1 

£ 

1 


1 


- 0. (35) 


1 

£ 


As one can see, we get a result which is typical for 
the perturbation process: the fourfold degenerate level e — 0 
splits partially and the non-degenerate level moves. At the 
same time level £ = 0 still remains threefold degenerate. That 
is why there is no arising gap between the valence and con¬ 
duction bands. The level, which is split from the originally 
fourfold degenerated level e = 0 is the impurity level in the 
conduction band, that is the resonance level. 


For other values of x the equations have the similar form. 
In all cases the equations of the sixth degree are to be reduced 
to equations of the third degree and are solved analytically. 
The results of these calculations are: 

x = ^ : £1,2 = 0; £3 = A 

£4,5,6 are determined by formula (33), but with other c p, 
namely 
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X = 2’ 


cos ip 

(A 2 

(SiC) 



£1,2 — 


£3,4 — 


£ 5,6 = 


(36) 


(37) 


A ± VA 2 + 36 


x = - : ei = 0; 

£ 2,3 = A; 


^ 4 , 5,6 — 


A - \/A 2 -j-27 cos ( ^ + 


n = 0,1,2... 

where <p is the same as in (36) 


(38) 


5 

X = 6 ' 


£1,2,3 = A; 
2 

^ 4 , 5,6 - 3 


A + \/A 2 + 27 cos ( ^ + n 


where 


A(A 2 - |3 4 


COS (Z> — 

(A 2 + 27)^ 
n = 0,1, 2 ... 

(39) 

The dependence of the energy bands of Ci- X Si x on x is 
shown in form of a diagram in Fig. 14. This diagram is ob¬ 
tained from formulas (33),(34),(36)-(39) for A = 3.5 which 
corresponds roughly to the value of A for Si in graphene. 

The main feature of this band structure rises from forbid¬ 
den band with the typical tunnel form. From this one can 
suggest, that in principle it is possible to create a tunnel diode 
on the basis on a Ci_ x Si x layer. In the ground of such a diode 
can be a tunnelling between Si- and C-domains of the layer 
in which the concentration of x changes with the coordinate 
x. That is, if in Fig 14 x is not only the concentration, but as 
well the coordinate x. 

Next we study the clusters of type II in which impurity 
atoms interact not only with basic atoms of the layer but also 
with each other. The latter is possible, if the impurities are 
in the lattice points of different sublattice: A and B, because 
the atoms of the same sublattice do not interact with each 
other in the nearest neighbors approximation. 

In case of x = | the cluster of type II can be obtained from 
the cluster in Fig. 13(a), if place the atoms Si in the lattice 
points, for example 2 and 5 instead of 1 nad 3. C-atoms with 
x = | in the cluster in Fig. 13(c) can be treated the same way. 
As a result we obtain following solutions for clusters of type 
II: 


x = ^ : £ 1,2 - 0; 

£ 3,4 =d (a - 3 ± v/A 2 + 2A + 9 ) ; (40) 

£5,6 = ^ (A + 3 ± VA 2 - 2A + 9) . 



FIG. 14: Energy spectrum of the mono-atomic layer 
Ci- X Si x . checkered: valence band; striped: conduction 
band; plain: forbidden band 


x = 3 : £ 1,2 = A; 

£3,4 = 1 (a - 3 ± VA 2 - 2A + 9 ) ; (41) 

£5,6 — ^ + 3 ± A 2 2 A 9^ . 

From formulas 40 and 41 one can obtain the energy of im¬ 
purity states, which at A = 3.5 can be found most within 
valance and conduction bands. The single level which can be 
found in the forbidden band appears when x = |. This level 
can take part in tunnelling, which we discussed above. 


IX. CONCLUSIONS 

In this paper, the CC approach has been applied to some 
problems of calculation of the grpahene energy spectrum. 
This approach seems to be especially appropriate both for 
problems of the breach of the graphene’s periodical structure 
(point defects, boundaries) and for problems besides mono- 
layer graphene (bilayer graphene, compound Ci_ x Si x ). 

As for the precision of the approach, it may easily be in¬ 
creased by considering clusters of larger sizes and interaction 
not only with the nearest neighbors. Besides, it is not diffi¬ 
cult to study the cluster with a distortion of the lattice near 
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the point defect. One can also consider the conglomerates of 
point defects in a graphene lattice. 

Except various problems with the calculation of electron en¬ 
ergy spectrum, the closed cluster approach may as well be ap¬ 
plied to calculate of the vibration spectrum in graphene, like¬ 
wise as it is done for one-dimensional and three-dimensional 
crystals in [18], [19]. 
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